Figure 12.2

ω=2 π

(Subscript[ω, 0])=1.5 ω

β=(Subscript[ω, 0])/4

γ=0.2

Fig12p2=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==0,ϕ'[0]==0},ϕ,{t,0,6}]

Plot[Evaluate[ϕ[t]/.Fig12p2],{t,0,6},Ticks->{{1,2,3,4,5,6},{-0.3,0.3}},PlotRange->All,AxesStyle->Thick,PlotStyle->Thick,LabelStyle->{Bold,Medium}]

Figure 12.3

γ=0.9

Fig12p3a=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==0,ϕ'[0]==0},ϕ,{t,0,6}]

Plot[Evaluate[ϕ[t]/.Fig12p3a],{t,0,6},PlotRange->All,Ticks->{{1,2,3,4,5,6},{-π/2,π/2}},AxesStyle->Thick,PlotStyle->Thick,LabelStyle->{Bold,Medium}]

Plot[{Evaluate[ϕ[t]/.Fig12p3a],2.7Cos[ω t-1.76]},{t,4.9,6.1},Ticks->{{5,6},{-π/2,π/2}},PlotRange->All, PlotStyle->{Thick,Dashed},AxesStyle->Thick,LabelStyle->{Bold,Large},AxesOrigin->{4.9,0}]

Plot[{Evaluate[ϕ[t]/.Fig12p3a],2.7Cos[ω t-1.76],-0.2Cos[3*(ω t-1.76)]},{t,4.9,6.1},Ticks->{{5,6},{-π/2,π/2}},PlotRange->All, PlotStyle->{Thick,Dashed},AxesStyle->Thick,LabelStyle->{Bold,Large},AxesOrigin->{4.9,0}]

Figure 12.4

γ=1.06

Fig12p4=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==0,ϕ'[0]==0},ϕ,{t,0,15}]

Plot[Evaluate[ϕ[t]/.Fig12p4],{t,0,15},PlotRange->All, Ticks->{{3,6,9,12,15},{0,2 Pi,4 Pi}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Figure 12.5

γ=1.073

Fig12p5a=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==0,ϕ'[0]==0},ϕ,{t,0,30}]

Plot[Evaluate[ϕ[t]/.Fig12p5a],{t,0,30},PlotRange->All, Ticks->{{5,10,15,20,25,30},{-2 π,0,2 Pi,4 Pi}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Plot[Evaluate[ϕ[t]/.Fig12p5a],{t,20,30},PlotRange->{-8.7,-7.8}, Ticks->{{20,22,24,26,28,30},{-8.5,-8.0}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Figure 12.6

γ=1.077

Fig12p6=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==0,ϕ'[0]==0},ϕ,{t,0,15}]

Plot[Evaluate[ϕ[t]/.Fig12p6],{t,0,15},PlotRange->All, Ticks->{{2,4,6,8,10,12,14},{0,2 Pi,4 Pi}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Figure 12.7

γ=1.077

Fig12p7=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==0,ϕ'[0]==0},ϕ,{t,0,15}]

Fig12p7p2=NDSolve[{ϕ1''[t]+2 β ϕ1'[t]+(Subscript[ω, 0])^2 Sin[ϕ1[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ1[0]==-π/2,ϕ1'[0]==0},ϕ1,{t,0,15}]

Plot[{Evaluate[ϕ[t]/.Fig12p7],Evaluate[ϕ1[t]/.Fig12p7p2]},{t,0,15},PlotRange->All, Ticks->{{2,4,6,8,10,12,14},{0,2 Pi,4 Pi}},PlotStyle->{Dashed,Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Figure 12.8

γ=1.06

Fig12p8a=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ,{t,0,40}]

Plot[Evaluate[ϕ[t]/.Fig12p8a],{t,0,10},PlotRange->All, Ticks->{{2,4,6,8,10},{-2,0,2}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Plot[Evaluate[ϕ[t]/.Fig12p8a],{t,28,40},PlotRange->{-1.7,-2.5}, Ticks->{{28,30,32,34,36,38,40},{-2.5,-2}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

γ=1.078

Fig12p8a=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ,{t,0,40}]

Plot[Evaluate[ϕ[t]/.Fig12p8a],{t,0,10},PlotRange->All, Ticks->{{2,4,6,8,10},{-2,0,2}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Plot[Evaluate[ϕ[t]/.Fig12p8a],{t,28,40},PlotRange->{-1.7,-2.5}, Ticks->{{28,30,32,34,36,38,40},{-2.5,-2}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

γ=1.081

Fig12p8a=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ,{t,0,40}]

Plot[Evaluate[ϕ[t]/.Fig12p8a],{t,0,10},PlotRange->All, Ticks->{{2,4,6,8,10},{-2,0,2}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Plot[Evaluate[ϕ[t]/.Fig12p8a],{t,28,40},PlotRange->{-1.7,-2.5}, Ticks->{{28,30,32,34,36,38,40},{-2.5,-2}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

γ=1.0826

Fig12p8a=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ,{t,0,40}]

Plot[Evaluate[ϕ[t]/.Fig12p8a],{t,0,10},PlotRange->All, Ticks->{{2,4,6,8,10},{-2,0,2}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Plot[Evaluate[ϕ[t]/.Fig12p8a],{t,28,40},PlotRange->{-1.7,-2.5}, Ticks->{{28,30,32,34,36,38,40},{-2.5,-2}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Figure 12.10

γ=1.105

Fig12p10=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ,{t,0,30}]

Plot[Evaluate[ϕ[t]/.Fig12p10],{t,0,30},PlotRange->All, Ticks->{{5,10,15,20,25,30},{-π,π}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Figure 12.11

γ=0.1

Fig12p11a=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ,{t,0,10},AccuracyGoal->100,WorkingPrecision->30]

Fig12p11b=NDSolve[{ϕ1''[t]+2 β ϕ1'[t]+(Subscript[ω, 0])^2 Sin[ϕ1[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ1[0]==-π/2-0.1,ϕ1'[0]==0},ϕ1,{t,0,10},AccuracyGoal->100,WorkingPrecision->30]

Plot[Log[10,Abs[Evaluate[ϕ[t]/.Fig12p11a]-Evaluate[ϕ1[t]/.Fig12p11b]]],{t,0,10},PlotRange->{-12,0}, Ticks->{{2,4,6,8,10},{-12,-10,-8,-6,-4,-2,0}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Figure 12.12

γ=1.07

Fig12p12a=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ,{t,0,40},AccuracyGoal->200,WorkingPrecision->18]

Fig12p12b=NDSolve[{ϕ1''[t]+2 β ϕ1'[t]+(Subscript[ω, 0])^2 Sin[ϕ1[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ1[0]==-π/2-0.1,ϕ1'[0]==0},ϕ1,{t,0,40},AccuracyGoal->200,WorkingPrecision->18]

Plot[Log[10,Abs[Evaluate[ϕ[t]/.Fig12p12a]-Evaluate[ϕ1[t]/.Fig12p12b]]],{t,0,40},PlotRange->{-8,0}, Ticks->{{5,10,15,20,25,30,35,40},{-8,-6,-4,-2,0}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Figure 12.13

γ=1.105

Fig12p13a=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ,{t,0,17},AccuracyGoal->200,WorkingPrecision->25]

Fig12p13b=NDSolve[{ϕ1''[t]+2 β ϕ1'[t]+(Subscript[ω, 0])^2 Sin[ϕ1[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ1[0]==-π/2-10^(-4),ϕ1'[0]==0},ϕ1,{t,0,17},AccuracyGoal->200,WorkingPrecision->25]

Plot[Log[10,Abs[Evaluate[ϕ[t]/.Fig12p13a]-Evaluate[ϕ1[t]/.Fig12p13b]]],{t,0,17},PlotRange->{-6,1}, Ticks->{{2,4,6,8,10,12,14,16},{-6,-5,-4,-3,-2,-1,0,1}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Figure 12.14

γ=1.13

Fig12p14a=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ,{t,0,20},AccuracyGoal->200,WorkingPrecision->25]

Plot[Evaluate[ϕ[t]/.Fig12p14a],{t,0,15},PlotRange->All, Ticks->{{5,10,15},{0,2 π}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Fig12p14b=NDSolve[{ϕ1''[t]+2 β ϕ1'[t]+(Subscript[ω, 0])^2 Sin[ϕ1[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ1[0]==-π/2-10^(-3),ϕ1'[0]==0},ϕ1,{t,0,20},AccuracyGoal->200,WorkingPrecision->25]

Plot[Log[10,Abs[Evaluate[ϕ[t]/.Fig12p14a]-Evaluate[ϕ1[t]/.Fig12p14b]]],{t,0,20},PlotRange->{-9,0}, Ticks->{{0,5,10,15,20},{-9,-6,-3,0}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Figure 12.15

γ=1.503

Fig12p15a=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ,{t,0,25},AccuracyGoal->200,WorkingPrecision->23]

Plot[Evaluate[ϕ[t]/.Fig12p15a],{t,0,25},PlotRange->All, Ticks->{{5,10,15,20,25},{-20 π,-10 π,0}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Fig12p15b=NDSolve[{ϕ1''[t]+2 β ϕ1'[t]+(Subscript[ω, 0])^2 Sin[ϕ1[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ1[0]==-π/2-10^(-3),ϕ1'[0]==0},ϕ1,{t,0,25},AccuracyGoal->200,WorkingPrecision->23]

Plot[Log[10,Abs[Evaluate[ϕ[t]/.Fig12p15a]-Evaluate[ϕ1[t]/.Fig12p15b]]],{t,0,15},PlotRange->{-4.5,1.75}, Ticks->{{0,5,10,15},{-4,-3,-2,-1,0,1}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Figure 12.16

γ=1.503

Fig12p16a=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ,{t,0,25},AccuracyGoal->200,WorkingPrecision->23]

Fig12p16b=NDSolve[{ϕ1''[t]+2 β ϕ1'[t]+(Subscript[ω, 0])^2 Sin[ϕ1[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ1[0]==-π/2-10^(-3),ϕ1'[0]==0},ϕ1,{t,0,25},AccuracyGoal->200,WorkingPrecision->23]

Plot[{Evaluate[ϕ[t]/.Fig12p16a],Evaluate[ϕ1[t]/.Fig12p16b]},{t,0,25},PlotRange->All, Ticks->{{5,10,15,20,25},{-20 π,-10 π,0}},PlotStyle->{Thick},AxesStyle->Thick,LabelStyle->{Bold,Large}]

Figure 12.19

γ=1.4

Fig12p19a=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ,{t,0,10},AccuracyGoal->200,WorkingPrecision->30]

Plot[Evaluate[ϕ[t]/.Fig12p19a],{t,0,10},PlotRange->All, Ticks->{{5,10},{-20 π,-10 π,0}}]

Fig12p19b=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ',{t,0,10},AccuracyGoal->200,WorkingPrecision->30]

Plot[Evaluate[ϕ'[t]/.Fig12p19b],{t,0,10},PlotRange->All, Ticks->{{5,10},{-20 ,-10 ,0,10}}]

Figure 12.20

γ=0.6

Fig12p20=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ,{t,0,6},AccuracyGoal->200,WorkingPrecision->30]

Plot[Evaluate[ϕ[t]/.Fig12p20],{t,0,6},PlotRange->All, Ticks->{{1,2,3,4,5,6},{-π/2,π/4}}]

Figure 12.21

γ=0.6

Fig12p21a=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ,{t,0,20},AccuracyGoal->200,WorkingPrecision->25]

Fig12p21b=NDSolve[{ϕ1''[t]+2 β ϕ1'[t]+(Subscript[ω, 0])^2 Sin[ϕ1[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ1[0]==-π/2,ϕ1'[0]==0},ϕ1',{t,0,20},AccuracyGoal->200,WorkingPrecision->25]

Flatten[{Evaluate[ϕ[2]/.Fig12p21a],Evaluate[ϕ1'[2]/.Fig12p21b]}]

ParametricPlot[Flatten[{Evaluate[ϕ[t]/.Fig12p21a],Evaluate[ϕ1'[t]/.Fig12p21b]}],{t,0,20},PlotRange->All, Ticks->{{-π/2,-π/4,π/4},{-4,4}},AspectRatio->1/GoldenRatio]

ParametricPlot[Flatten[{Evaluate[ϕ[t]/.Fig12p21a],Evaluate[ϕ1'[t]/.Fig12p21b]}],{t,5,20},PlotRange->All, Ticks->{{-π/2,-π/4,π/4},{-4,4}},AspectRatio->1/GoldenRatio]

Figure 12.22

γ=0.6

Fig12p22a=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==0,ϕ'[0]==0},ϕ,{t,0,20},AccuracyGoal->200,WorkingPrecision->25]

Fig12p22b=NDSolve[{ϕ1''[t]+2 β ϕ1'[t]+(Subscript[ω, 0])^2 Sin[ϕ1[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ1[0]==0,ϕ1'[0]==0},ϕ1',{t,0,20},AccuracyGoal->200,WorkingPrecision->25]

ParametricPlot[Flatten[{Evaluate[ϕ[t]/.Fig12p22a],Evaluate[ϕ1'[t]/.Fig12p22b]}],{t,0,20},PlotRange->All, Ticks->{{-π/2,-π/4,π/4},{-4,4}},AspectRatio->1/GoldenRatio]

ParametricPlot[Flatten[{Evaluate[ϕ[t]/.Fig12p22a],Evaluate[ϕ1'[t]/.Fig12p22b]}],{t,5,20},PlotRange->All, Ticks->{{-π/2,-π/4,π/4},{-4,4}},AspectRatio->1/GoldenRatio]

Figure 12.23 Period 2

γ=1.078

Fig12p23a=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ,{t,0,60},AccuracyGoal->200,WorkingPrecision->14]

Fig12p23b=NDSolve[{ϕ1''[t]+2 β ϕ1'[t]+(Subscript[ω, 0])^2 Sin[ϕ1[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ1[0]==-π/2,ϕ1'[0]==0},ϕ1',{t,0,60},AccuracyGoal->200,WorkingPrecision->14]

ParametricPlot[Flatten[{Evaluate[ϕ[t]/.Fig12p23a],Evaluate[ϕ1'[t]/.Fig12p23b]}],{t,20,60},PlotRange->All, Ticks->{{-π/2,π/2},{-10,10}},AspectRatio->1/GoldenRatio]

Figure 12.23 Period 4

γ=1.081

Fig12p23a4=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ,{t,0,60},AccuracyGoal->200,WorkingPrecision->14]

Fig12p23b4=NDSolve[{ϕ1''[t]+2 β ϕ1'[t]+(Subscript[ω, 0])^2 Sin[ϕ1[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ1[0]==-π/2,ϕ1'[0]==0},ϕ1',{t,0,60},AccuracyGoal->200,WorkingPrecision->14]

ParametricPlot[Flatten[{Evaluate[ϕ[t]/.Fig12p23a4],Evaluate[ϕ1'[t]/.Fig12p23b4]}],{t,20,60},PlotRange->All, Ticks->{{-π/2,π/2},{-10,10}},AspectRatio->1/GoldenRatio]

Figure 12.24

γ=1.105

Fig12p24a=NDSolve[{ϕ''[t]+2 β ϕ'[t]+(Subscript[ω, 0])^2 Sin[ϕ[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ[0]==-π/2,ϕ'[0]==0},ϕ,{t,0,200},AccuracyGoal->200,WorkingPrecision->9]

Fig12p24b=NDSolve[{ϕ1''[t]+2 β ϕ1'[t]+(Subscript[ω, 0])^2 Sin[ϕ1[t]]==γ (Subscript[ω, 0])^2 Cos[ω t],ϕ1[0]==-π/2,ϕ1'[0]==0},ϕ1',{t,0,200},AccuracyGoal->200,WorkingPrecision->9]

ParametricPlot[Flatten[{Evaluate[ϕ[t]/.Fig12p24a],Evaluate[ϕ1'[t]/.Fig12p24b]}],{t,14,21},PlotRange->All, Ticks->{{-1,1},{-10,10}},AspectRatio->1/GoldenRatio]

ParametricPlot[Flatten[{Evaluate[ϕ[t]/.Fig12p24a],Evaluate[ϕ1'[t]/.Fig12p24b]}],{t,14,94},PlotRange->All, Ticks->{{-1,1},{-10,10}},AspectRatio->1/GoldenRatio]